Longitudinal and transverse velocity scaling exponents from merging of the Vortex 

filament and Multifractal models. 
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We suggest a simple explanation of the difference between transverse and longitudinal scaling 
exponents observed in experiments and simulations. Based on the Vortex filament model and 
Multifractal conjecture, we calculate both scaling exponents for any n without any fitting parameters 
and ESS anzatz. The results are in very good agreement with the data of simulations. 



Numerous direct numerical simulations (DNS) of hy- 
drodynamical turbulent flow, as well as high-resolution 
experiments, have been performed during last twenty 
years 1]. The positive result of these investigations 
is that numerical and experimental approaches demon- 
strate a very good agreement. In particular, velocity 
statistics calculated from DNS practically coincides with 
that obtained from experiments [l[ . This justifies the ac- 
curacy of both kinds of results. But there still is a lack of 
physical understanding of the processes occuring in a tur- 
bulent flow and contributing to statistics. Which of the 
infinite number of solutions to the NSE are responsible 
for the observed intermittent behavior, and why? 

Okamoto et al. Q claimed that the most part of the 
flow's helicity and, in some sense, the most part of infor- 
mation about the flow are contained in the same small 
part of the liquid volume spread more or less uniformly 
throughout the liquid. These structures are stable, their 
lifetime exceeds many times the largest eddy turnover 
time in the flow. 

This seems us to be a very important step. Proba- 
bly, these elongated regions of high vorticity are just the 
structures that determine the small-scale statistics of a 
turbulent flow. Thus, it would be useful to elaborate 
some physical understanding of these objects. 

Such an attempt was performed in The main 

idea is that in the regions where vorticity is high (vor- 
tex filaments), the vorticity itself stabilizes the motion. 
So, despite the stochastic large-scale forces that give en- 
ergy to a filament, the motion inside the filament has 
an essential non-stochastic component. This may be the 
reason of stability of filaments. The random large-scale 
forces then act, on average, as a stretching force. Stretch- 
ing (not breaking!) of vortices is the mechanism that 
provides the observed energy flux from larger to smaller 
scales, and observed statistical properties of turbulence. 
This is the main difference of our approach from the Kol- 
mogorov's approach developed in the K41 theory [111 ]. 

The relation between longitudinal and transversal Eu- 
ler velocity structure functions is one of the problems 
that have had no answer up to now. These functions are 
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Here Av = v(r + 1) — v(r)is velocity difference between 
two near points separated by 1, and the average is taken 
over all pairs of points separated by given I. 

Experiments and DNS show that inside the inertial 
range the structure functions obey scaling laws, 

#(0 cx # , srU(i) oc A . 

The scaling exponents £n> Crl are believed to be inde- 
pendent on conditions of the experiment. They are in- 
termittent, i.e., Cn/ n an d Cn/n are decreasing functions 
of n. 

The question whether the two scaling exponents co- 
incide in isotropic and homogeneous turbulence or not, 
is open. There is an exact statement that cl = Cn f° r 
2 and 3 ll|. On the other hand, modern exper- 
and numerical simulations show 
some significant difference between y| and at higher 
n. But the proponents of the equality argue that the dif- 
ference may result from finite Reynolds number effects 
or anisotropy [Hj]. 
The aim of the paper is to analyze the question by 
means of the Vortex filament model. We show that the 
divergence of the scaling exponents can be understood 
and calculated from general physical considerations using 
the tools provided by the Vortex filament model and the 
well-known Multifractal (MF) conjecture. 

Before we consider the problem from the viewpoint of 
the Vortex filament model, we note that isotropic and ho- 
mogeneous medium does not necessarily imply = £„. 
For a simple illustration, consider a gas of tops, each ro- 
tating around its own axis, the axes directions distributed 
randomly. Although locally a strong asymmetry could be 
found near each top, the whole picture remains isotropic. 
The Vortex filament model provides analogous situation 
that might be in real turbulence. 

In [l0( we suggested to join our model of vortex fil- 
aments with the MF model to get velocity scaling ex- 
ponents. We now remind briefly some ideas of the two 
models and the results of [lfj. Then we proceed to the 
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difference between the longitudinai and transverse scal- 
ing exponents. 

The MF model implies that the determinative contri- 
bution to velocity structure functions is given by the so- 
lutions (regions) where 

Av{l) = |v(r + l) -v(r)| ~ 

So, to calculate structure functions, it is enough to con- 
sider only a set of scaling solutions that can be numbered 
by h. This property is called local scale invariance. 

The Large Fluctuations Theorem states that the prob- 
ability of measured velocity difference Av(l) to have the 
scaling h is a power-law function of I: 

P oc l 3 ~ D ^ 

Knowing D(h), one could in principle calculate all struc- 
ture functions: 

(Av n ) = [ 5v n Pdh = I r h l 3 ~ D ^d^(h) 



Here du(h) is the measure that is responsible for relative 
weights of different values h. In the limit I — >■ 0, only 
the smallest exponent contributes to the integral; it then 
follows 



min (nh 

h 



■3-D(h)) 
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Without loss of generality, one can treat D(h) as a con- 
cave function. 

As follows from |T]), the point ho where D(h) reaches its 
maximum corresponds to n = 0: the requirement Co = 
leads to D(ho) = 3. 

The definitional domain of D is restricted by some 
hmin and h max . Since we are interested in positive val- 
ues n, we will hereafter restrict ourselves by h < ho- 
What about h m in> from the condition ll| < an d 
the finiteness of the Mach number (or from the absence 
of singularities of Av as I —> 0) it follows (' n > for any 
n. The minimum in ([I]) is reached at h n = Q' n . Since 
negative values of the derivative are forbidden, we get 
hmin > 0. The behavior of D(h) near h m i n determines 
the behavior of at 

n ~ 7i* = D'(h min ) . 

For larger n, the minimum in ([1]) is reached at the bound- 
ary h = h m i n . Thus, the behavior of £ n at very large n 
depends only on the value of h m i n : the asymptotic be- 
havior of C« for all non-zero h m i n is 



3-D(h 
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0, the possible growth of £ n at n > n* is steeper 



If h m 

than linear; if D{h m i n ) is finite, we get Q n 



We now proceed to the Vortex filament model. In Q 
we derived the equation 



Pij 



(3) 



Here u> is the vorticity of the flow, uj = V x v; . 
is the pressure hessian. This equation is the direct con- 
sequence of the inviscid limit of the NSE and describes 
the evolution of vorticity along the trajectory of a liquid 
particle. It was for the first time obtained in 22J. The 
main assumption of the theory is that inside vortex fil- 
aments the 'longitudinal' (in relation to us) part of the 
hessian doesn't depend on the local vorticity, and is de- 
termined by the large-scale component of a flow. 1 This 
makes © a linear stochastic equation in relation to us. In 
3] we wrote the corresponding equation for probability 
density function (PDF), f(uj,uj,t). Solving the equation 
and integrating over all angles and over w, we found the 
intermediate asymptotic solution for the PDF at uj — > oo: 

P{uj) oc 1/uj 4 

In [l(| we used it to derive the condition 



min(3/i + 2 - D(h)) = 



(4) 



for the function D(h). As follows from (TTJ), this condition 
is just equivalent to 

C3 = 1 (5) 
which is the well-known Kolmogorov's '4/5 law'. 

The MF model is a dimensional theory, so it does 
not distinguish longitudinal and transverse scaling ex- 
ponents. Also the condition (fj| is valid for both of them. 
To describe the difference, it is natural to use two differ- 
ent functions D{h) for and CrJ 



17] 



But what is the 
relation between them, why does the difference appear? 

In 0, i we introduced a simple model of an axially 
symmetric vortex filament: 

v r = a(t)r , v z = b(t)z , = w(t)r (6) 

From the incompressibility condition and the Euler equa- 
tion it then follows 

2a + b = , a + a 2 - uj 2 = -P a (t) , 

w + 2acj = 0, b + b 2 = -P 2 {t), (7) 

p= \P x {t)r 2 + \P 2 {t)z 2 

Taking the second derivative of u), we get 

w = -P 2 (t)w , 



ln>n* 



const. 



[1] An indirect confirmation of this statement is presented in |2.'J . 
where the maximal eigenvalue of the pressure hessian is shown 
to be orthogonal to vorticity. 
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which corresponds to (|3]) . In accordance with the as- 
sumption discussed above, let P?{t) be a random func- 
tion independent of u>. Then all moments of u>(t) increase 
exponentially as a function of time, a ~ di/u) <C u> for 
large t. 

It is easy to calculate that in this model S„ oc (u> n )l n , 
Sn oc {a n )l n . In jg| we discussed the connection between 
t and I and showed that t oc — InZ. Hence, in this model 

S^ 3> Sn, Cit C«- 

The equality of scaling exponents can be restored if we 
take the other branch of the same filament into account, 
assuming that the filament is closed and 'turns back' at 
very large \z\. Then there are two filaments with equal 
and opposite vorticities. The second branch produces 
a perturbation of velocity Sv ~ aw, where the coeffi- 
cient a is inversely proportional to the distance between 
the two branches. (This is just analogous to an axial 
electric current producing magnetic field in vacuum; one 
can write the corresponding solution to the Euler equa- 
tion by developing the perturbation as series in a.) The 
perturbation violates the axial symmetry, so the scal- 
ing exponents become equal. However, the pre-exponent 
of Sn is much smaller than that of S„ if ex. is small. 
The value of a increases as we approach the turnover 
point, and the two branches become nearer. This gives 
us a hint that the regions where a filament is strongly 
curved may make a small contribution to S x but con- 
tribute much to Sn. This effect is stronger for large n, 
since Sjl ~ ct n S^- So, for infinitely large n we can ex- 
pect that transverse structure functions are dominated 
by extremely stretched-out and very thin, roughly axi- 
ally symmetric filaments with very high vorticities, while 
for longitudinal structure functions one needs extremely 
curved parts of these filaments. 

The discussed example is one of numerous solutions 
to the Euler equation that are presented in a turbulent 
flow. The observed exponents and £ n are produced by 
contributions of many filaments. Different filaments (and 
different parts of them) make contributions to different n 
(or, in terms of MF theory, to different h). This picture 
is very complicated, but, thanks to the MF model, we 
don't need to know it in details. Approximating D(h) by 
a second-order polynomial [13] and making use of (j2J , we 
only have to fix one more point of the curve. 

For this purpose, we consider the minimal possible h = 
h m in- This corresponds to n — > oo. In [l0( we have 
shown that h m - m = 0: values h < are forbidden by the 
condition d£ n /dn > (see p. 2). 

The existence of ft. = is proved by the possibility of a 
cylindric filament with rotating velocity independent of r: 

v = [e„r/r] (8) 

One can check that in this extreme case, indeed, Sv(l) ~ 
1°. The velocity {SJ satisfies the Euler equation with 



pressure logarithmically divergent; this means that for 
any positive h, there exists a corresponding velocity dis- 
tribution with converging pressure. 

Calculation of the correlators directly from ((HJ gives 
under the limit n — > oo: 

on 

(|Au x l/i| n ) oc — I 2 , (\Avl/l\ n ) octi- 5/2 Z 2 (9) 
n 

The proportionality to I 2 is caused by the axial geometry 
of the filament (integrating rdr). This corresponds to the 
definition of D given in, e.g., [Ill ]: the probability that at 
least one of a pair of points would get inside the filament 
is proportional to I 2 . Thus, we get 

D x (0) = 1 (10) 

To find the difference between S„ and S n , we now remind 
that, at infinitely large n, S„ are contributed mostly by 
'cylindric' parts of filaments, while sj[ are dominated by 
point-like regions where filaments are bent very strongly. 
Actually, the behavior of the pre-exponents in (J9j) shows 
that the contribution of the 'extreme' filament to S x in- 
creases as n — y oo. To the opposite, its contribution to 
Sn for large n (i.e., small h) is very small. There must 
be other solutions to the Euler equation to determine 
the behavior of D"(h) for small h and, equivalently, to 
make the most contribution to sj[ for large n. These so- 
lutions correspond to 'strongly curved' extreme filament. 
To satisfy h — 0, velocity must be independent of r. It 
may, for example, take the form v = (v r (9),vg(8),0). 
Such a solution exists but it cannot be written analyt- 
ically. At 9 = 0, the radial velocity diverges weakly: 
vg oc — 9y/\n (1/9) , v r oc (1/9); the pressure diver- 
gence is p oc lnr, just as in (JSJ. Since 5v ~ 1° in the 
case, and averaging includes r 2 dr, the correlator is pro- 
portional to I 3 . This corresponds to 

D«(0)=0 (11) 

The difference between the boundary values (fTU)) and (fTTj) 
determines the difference between the functions D± and 
£>". 

We seek the solution D(h) in the simplest form 

D(h) =3-b(h- h Q ) 2 

and use ((3]) and (|10|) . (TTTT) to find the two unknown pa- 
rameters. We obtain two equations: 

^o-l(^) 2 -|' bU=3/(hlf for ^11 

Each of the equations has two roots. The bigger roots 
are non-physical, because in this case the curve £n would 
be constant already at n > — D' h (Q) = 2bho ~ 2, and 
(HI would not hold. For the second roots, in accordance 
with {T]), we have 

Cn = nho — n 2 /4b , n < n* = 26/iq (12) 
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FIG. 1: Euler longitudinal (upper branch) and transverse 
(lower branch) scaling exponents: the results of DNS from 
[ill (O) D) an d [H[ (A- v)> an d the prediction of the theory 
(p) (lines). 



Thus, we propose an explanation of the difference be- 
tween C| and based on the difference between the 
filaments that contribute to the two structure functions: 
roughly speaking, this is the difference between axially 
symmetric and strongly curved ones. This allows to find 
the values of ( n for very large n, and merging of the Vor- 
tex filament and Multifractal theories gives the whole 
functions. The obtained solutions ([L^ -lfTB" ]) fit very well 
the observed scaling exponents and C« , and we hope 
that they reveal the nature of the difference between lon- 
gitudinal and transverse structure functions. 

We thank Prof.A.V. Gurevich for his kind interest to 
our work. The work was partially supported by the RAS 
Program 'Fundamental Problems of Nonlinear Dynam- 
ics'. 



Substituting the values of ho,b, we get the scaling expo- 
nents: 



r - 

Sn 



0.367n- 1.12- 10 _ V , n < 16.3 ; 



3 , 



n > 16.3 



0.391n - 1.91 • 10~ 2 n 2 , n < 10.2 ; 



(13) 



2 , 



n > 10.2 



In Fig. 1 we compare this theoretical prediction with 
experimental data by [l6| and [l5| . We see that the result 
of our simple model is very close to the experimental 
results and lies inside the error bars of the experiments. 

Our simple model has one difficulty: the two parabo- 
las plotted in Fig. 1 coincide at n — and n — 3, hence 
they do not coincide at n = 2. This contradicts to the 
exact theoretical statement that £2 = X2- However, this 
difficulty is caused by postulating the simplest parabolic 
shape for D(h). It can easily be solved by adding one 
more parameter and assuming D(h) to be a cubic poly- 
nomial. Because of very small divergence (1.6 • 10~ 2 ) 
between £2 and %2 in Fig. 1, the coefficient by the eldest 
order would be very small (~ 10~ 4 ) . It would change very 
slightly (unnoticeable for an eye) the lines presented in 
Fig. 1. The only thing they may change significantly is 
the rate of approaching the constant at large (but still in- 
termediate) n (from 10 to 15, approximately). But this 
range of n is, anyway, badly described by the lowest- 
order polynomials: adding more degrees of freedom with 
very small coefficients, though unimportant for smaller 
n, would change the solutions for these n. However, the 
changes cannot be very big, since the exponents are still 
restricted by the values 2 and 3, respectively. 

One more comment is that, knowind D{h) 1 one can use 
the MF model to calculate, e.g., the PDF of velocity gra- 
dients or accelerations. In [l[ it is shown that once D(h) 
fits Cn well, it would also fit well the other quantities. 
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